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A comprehensive study of the electronic states at the 4s+5s asymptote in KRb is presented. 
Abundant spectroscopic data on the a'^E"*" state were coUected by Fourier-transform spectroscopy 
which allow to determine an accurate experimental potential energy curve up to 14.8 A . The existing 
data set (C. Amiot et al. J. Chem. Phys. 112, 7068 (2000)) on the ground state X^E^ was extended 
by several additional levels lying close to the atomic asymptote. In a coupled channels fitting routine 
complete molecular potentials for both electronic states were fitted. Along with the line frequencies 
of the molecular transitions, recently published positions of Feshbach resonances in ""^K and **''Rb 
mixtures (F. Ferlaino et al. Phys. Rev. A 74, 039903 (2006)) were included in the fit. This makes 
the derived potential curves capable for an accurate description of observed cold collision features 
so far. Predictions of scattering lengths and Feshbach resonances in other isotopic combinations are 
reported. 

PACS numbers: 31.50.Bc, 33.20.Kf, 33.20.Vq, 33.50.Dq 



I. INTRODUCTION 

Along with NaK and NaRb, the KRb molecule is one 
of the most extensively studied among the mixed alka- 
lis. The ground state X^E"*", the lowest singlet states and 
some of the triplet states were studied by various spec- 
troscopic techniques. It is therefore surprising that the 
triplet state a^S"*" was not characterized experimentally 
so far. 

In recent years, understanding interactions between 
different alkali-metal species in their singlet and triplet 
ground states of the pairs became important due to the 
experiments in mixed atomic traps. This stimulated in- 
vestigations of the molecular states as well. Using molec- 
ular spectroscopy, several ground state asymptotes were 
examined: NaRb 3, NaCs f?], LiCs 0]. In the present 
study we applied the same technique for KRb. This 
molecule is particularly interesting since several groups 
are operating two-species traps, cold KRb molecules have 
been already formed, positions of several Feshbach res- 
onances in ^"K®''Rb have been measured [1, 0, 0] and 
new ones are still being searched for. The study of 
the (4s)K-|-(5s)Rb asymptote has therefore several strong 
motivations: (i) to obtain complete and accurate poten- 
tials for the X^I]+ and a^^S^ states, to be able to model 
cold collisions; (ii) to prove of the validity of the iso- 
topic scaling rules, i.e. applying the potentials derived 
for the main isotopic combination '^^K^'^Rb to reproduce 
the experimental observations (including the Feshbach 
resonances) in other isotopes of KRb; (iii) to predict new 
Feshbach resonances that can be later verified and pre- 
cisely measured experimentally with much reduced effort 
for the search procedure. 



Along with the theoretical calculations of the long 
range dispersion coefficients 0, H, H, there is only 
one theoretical work yet devoted to the study of the elec- 
tronic structure of KRb [ll[ . 

The ground state X-^E+of KRb was first studied for 
v" < 44 at high resolution by Ross and coworkers ■ A 
significant step towards the asymptote was done in Ref. 
[l3|, where an accurate potential energy curve for v" up 
to 87 was derived. In both investigations the laser in- 
duced fluorescence to the KRb ground state was recorded 
by exciting the A^S"*" state by radiation of a Ti: S app hire 
laser. The B{iyU state was studied in Refs. [13, El Ell- 
Numerous perturbations were observed due to mixing 
with the nearby lying triplet states. Therefore in or- 
der to study the a'^E"'" state it was expected to follow 
straightforward the strategy which was successfully ap- 
plied fot the investigation of the same states in NaRb [l| 
and NaCs [5], i.e. to use the perturbations in the B(l)^n 
state as a window to the triplet manifold. 



II. EXPERIMENT 

The KRb molecules were produced in a stainless steel 
heat-pipe similar to that developed for NaRb [l| and 
NaCs 12;]. Metallic K and Rb were mixed in an ap- 
proximate mass ratio 1:1. The heat-pipe was operated 
typically at a temperature of 550 K with 2-5 mbar of 
Ar as buffer gas. The Laser Induced Fluorescence (LIF) 
was excited with narrow band single mode lasers and 
resolved with a Fourier Transform Spectrometer (FTS) 
Bruker 120 HR. 

For excitation of the B (1)^11 state we used a Coher- 



2 



,'"=23 19 1615 14 12 11 10 9 



7 6 5 4 




10950 11000 

wave number, cm ' 



FIG. 1: (color online) A fragment of a triplet progression 
excited at 15059.018 cm~^ and originating from (ii', J' — 36) 
in the B - X system. The doublets of the progression are 
formed by transitions to levels with A^" = 36 and 38 of the 
a^E^ state. For strong lines (e.g. v"=6) coUisional satellites 
are also visible. 



ent 599 linear laser with DCM dye delivering a power 
of about 70 mW. In the region between 15000 cm~^ and 
15300 cm^-'^ we found a lot of transitions in the B-X band 
system. Along with them, transitions in the similar sys- 
tem of K2 were often seen with comparable intensity. 
Going to higher frequencies the transitions in K2 domi- 
nated those in KRb. It is worth noting here, that when 
we tested a heat-pipe filled only with potassium, the flu- 
orescence due to K2 at similar experimental conditions 
(temperature, laser frequency, laser power) was at least 
a factor of 3 stronger than in the KRb heat-pipe. Obvi- 
ously in case of mixtures the concentration of molecules 
cannot be derived simply from the saturated vapour pres- 
sure of the correspondent pure metals as often stated 
in spectroscopic papers. This is confirmed also by the 
fact that Rb2 fluorescence was not observed in the stud- 
ied spectral region, although Rb2 was reported to absorb 
there [l3,[il[l3. 

The fluorescence from the B state led usually back to 
the ground X^E"*" state. However when exciting a per- 
turbed level with triplet admixture, fluorescence down to 
the triplet state a'^S+ was seen as well. It appeared in 
the region of 11000-12000 cm~^ well separated from the 
B-X fluorescence. An example is given in Fig. [Uwhere 
the fluorescence progression extends close to the dissoci- 
ation asymptote, which can be concluded from the large 
decrease of the vibrational spacing below 10950 cm~^. 
When resolved by the FTS with a resolution of typi- 
cally 0.03 cm^^ most lines to the a'^E^ state show the 
expected hyperflne structure (HPS) similar to that ob- 
served in Ref. P, • In those cases it could be explained 
by the Fermi contact interaction (20| using the atomic 
HFS parameters of Rb and K [2l|. However, contrary to 



the NaRb and NaCs cases P, , there was also a large 
number of progressions all lines of which show unusual 
HFS (see Fig. This we shall discuss in the following 
section IIIII 

The ground state X^E+ was studied in Ref. [l^l up 
to v"=87. The classical outer turning point of the last 
observed level is about 11.3 A. This value differs from 
that given in the title of Ref. [l^, namely 10 A, since 
the latter corresponds to the turning point of the level 
{v"=87, J"=0), but in fact the authors observed that vi- 
brational level for J"=51 which shifts the classical turn- 
ing point to a larger internuclear distance. This dis- 
tance, however is still too small in order to study the 
coupling between the a'^E"'" and X^E+ states due to the 
HFS mixing. Therefore, we revisited the X^E"*" state 
asymptote in order to extend the experimental informa- 
tion closer to the dissociation limit. For this purpose 
we excited the A^E+ - X^E+system using a Ti:Sapphire 
laser (from Tekhnoscan) pumped by a 10 W cw frequency 
doubled Nd:YAG laser (Coherent Verdi 10) and repeated 
as first steps measurements from [l3| for internal calibra- 
tion check. Then we extended the experimental data set 
by recording progressions from the same upper vibra- 
tional level but with different J'. It is worth noting that 
some of the lines of the progression reported in Ref. [l3| 
which go to v^=87 turned out to be overlapped with 
stronger lines from different progressions. Exciting other 
upper vibrational levels that would give fluorescence to 
high lying ground state levels was not applicable due to 
the presence of strong K2 fluorescence in the same spec- 
tral region. 

We searched also for transitions to the a^E"*" state due 
to the so called long-range-changeover which provided 
precious information in the case of NaRb [l| and LiCs . 
By exciting a B state level close to the atomic asymptote 
(in those cases (n^S)-l-(m^P)) we actually approach the 
range where Hund's case (c) starts to be more appro- 
priate than Hund's case (a) (due to the large flne struc- 
ture splitting of the lowest p state in Rb and Cs). As 
a consequence, a decay from such excited levels to both 
X^E+ {Q = 0+) and a'^E+ (fi = 1 component) states is 
possible, which cannot be explained within the Hund's 
case (a) or (b) models either due to selection rules or 
due to vanishing Franck-Condon factors for the case of 
local perturbations. The case of KRb, however, turned 
out to be more similar to that of NaCs 0], where the 
long-range changeover appeared rather fragmentary and 
provided only few transitions to near asymptotic levels 
in the X^E+ and a'^E+ states. 



III. ANALYSIS 

Our previous experience [l|, 0, Q has shown, that 
the states a'^E"'" in the heteronuclear alkali-metal dimers 
can be modeled within the experimental resolution in 
the Born-Oppenheimer approximation by a single poten- 
tial curve while the hyperfine splitting is well described 
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FIG. 2: (color online) Examples of "usual" and "unusual" HFS within the progressions to the a'^E^ state. On the left side 
the progression starts from an upper state level with J'=65. The doublets are formed by transitions to levels with N" = 67 
and 65, except those with v" = 7 and 8 where a weak transition to A'^' = 63 also appears. On the right side a transition from 
(«', J' = 39) is shown. Both insets give examples of hyperfine structure, details see text. 
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FIG. 3: (color online) Two hyperfine groups of closely spaced 
transitions to a^E^ state levels are shown for three differ- 
ent frequencies of the exciting laser. These transitions be- 
long to different progressions, i.e. they originate from differ- 
ent excited levels {v',J' = 9 and v',J' — 12). It is clearly 
seen that while the structure of the left transition (to the 
{v" — 6, A*'" = 9) a'^S'*' state level) remains almost unchanged 
and only the intensity varies, the overall pattern of the right 
transition (to the {v" — 6,N" — 11) level) changes; for the 
blue and red detuned excitation frequencies, different compo- 
nents of the hyperfine triplet are missing. 

within our resolution with the atomic hyperfine parame- 
ters. This is possible if the hyperfine structure does not 
depend on the internuclear separation. Only close to the 
atomic asymptote it is necessary to take into account the 
mixing with the ground states X-'^S"'". 

The usual appearance of the hyperfine splitting of the 



a'^I]+ levels is shown in Fig. [2] (inset of left graph). The 
splitting into three components is due to the interaction 
of the electrons with the nuclear spin of the Rb atoms. 
Each of these components is further split due to the in- 
teraction with the K nuclear spin. This interaction is 
small, contrary to the case of Na in NaRb for example 
and is below our experimental resolution. 

In the case of KRb we observed for the first time sev- 
eral progressions where the HFS differs from the expected 
one (see Fig. (2). We made special efforts to study these 
"unusual" features within the possibilities of our exper- 
imental setup. We found out that within a given pro- 
gression the hyperfine splitting is always independent 
of the rotational and vibrational quantum numbers of 
the a^S"*" (both for the usual and the unusual cases). 
At the same time transitions to one given a'^S"'' level 
coming from different excited levels can show different 
structure. Therefore, we concluded that it should be 
the structure of the upper state levels which is respon- 
sible for the deviations from the expected HFS. This 
idea was confirmed by the observation that the unusual 
HFS changes its appearance when the excitation laser 
is detuned within the Doppler profile of the molecular 
transition. In Fig. [3] we show results of such an ex- 
periment. The laser excites simultaneously two transi- 
tions, one ((«" = 0, J" = 9)— >(d', J' = 9)) gives rise to a 
progression with usual HFS (TV" — 7, 9, 11 are involved, 
IH) and the second {{v" = 0, J" = ll)-^(w', J' = 12), 
the assignment of J' is not certain, which is unimpor- 
tant for the analysis of the ground states) - with unusual 
one {N" = 11, 13 are involved). We chose a portion of 
the fluorescence spectrum where two lines, each belong- 
ing to one of the progressions, appear closely together. 
Then we set the laser frequency within the Doppler pro- 
file of the transitions at the three position (center and ± 
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0.018 cni~-'^and observed the change of the HFS of these 
hnes. It is clearly seen in Fig. [3] that while the structure 
of the line with N"=9 remains almost unchanged and 
varies only in intensity, the overall pattern of the line 
with N"=ll changes. For the blue and red detuned exci- 
tation frequencies, different components of the hyperfine 
triplet are missing. Moreover, from the relative changes 
of the intensities of the lines from one laser frequency to 
another we can expect that the effective widths of the ex- 
citation transitions are different. The transition produc- 
ing the progression with unusual HFS could be broader, 
which is probably due to HFS of the upper state level, be- 
ing comparable with the Doppler width of the transition. 
This hypothesis seems to be reasonable if we compare it 
with the HFS of the {2)^T,+ state in NaRb observed by 
Matsubara et al. in Ref. [1^. So, presently, we under- 
stand our observations in the following way: The HFS of 
the excited complex of singlet and triplet states in KRb 
is in some cases larger than the Doppler width of the 
transitions. The laser with given frequency excites a se- 
lected set of hyperfine components which in turn does 
not give fluorescence to the full multiplet of the a^S^ 
levels. Therefore, the observed structure of the spectral 
lines reflects the HFS both of the lower and the upper 
electronic states. 

As a center frequency for a given hyperfine multiplet 
we chose the frequency of the central component of the 
HFS triplet (see the left part of Fig. [2]) which is normally 
also the strongest. For the potential fit we used these fre- 
quencies and included the term energy of the upper level 
as a free parameter, common for all frequencies forming 
one progression. This selection is possible, because in a 
given progression the appearance of the HFS turned out 
to be independent of the vibrational and rotational quan- 
tum numbers within the experimental resolution. In case 
of progression with unusual HFS, we took a characteris- 
tic feature of the hyper fine group of lines (usually the 
strongest peak). In this way we were able to reproduce 
correctly the vibrational and rotational spacings of the 
a'^S"'" within a progression. We did not make efforts to re- 
produce the unusual HFS of the transitions, because our 
present understanding of the problem indicates that we 
need to include also the HFS of the upper state, the study 
of which is outside the scope of this work and would ask 
for a different experimental setup to get sufficient data 
for a convincing analysis. 



A. The a=*E+ state 

Establishing the rotational numbering of the triplet 
progressions was straightforward because almost always 
they were accompanied by the corresponding progression 
to the ground X^S"*" state. For finding the proper vibra- 
tional numbering we first decided to rely on the possi- 
bility of scaling the molecular constants of one electronic 
state with the ratio of the reduced masses of different 
isotopomers. Experimentally, progressions belonging to 



the molecule with the less abundant isotope * Rb were 
easy to distinguish by the broader HFS compared to 
®^Rb. In our initial analysis we had several progressions 
in 39K85Rb and one in ^SR^^Rb. The range of N" [H 
was from 31 to 68 and v" spanning a range of 15 vibra- 
tional quanta. Several fits of Dunham-type coefficients 
with different vibrational assignment including both iso- 
topomers indicated the most probable one which we took 
for further analysis of the experimental data and poten- 
tial energy curve (PEC) construction. 

Once the initial vibrational numbering was suggested, 
we fitted a potential curve for the a^S"*" state using a 
point-wise representation f23| and used it to assign new 
vibrational progressions. Periodically, as the size of the 
experimental data set significantly increased, the PEC 
was refined. We noticed, however, that for larger N" 
there was a growing systematic disagreement between 
the observed and the calculated transition frequencies, 
which reached about 0.05 cm~^ for the highest N". Af- 
ter all experimental data were assigned, this problem re- 
mained. Similarly to the a^S"*" state in NaRb [Tl|, we 
can exclude the possibility for significant spin-rotational 
and second order spin-orbit effects in the a'^S^ state of 
KRb. Therefore, we searched for a connection between 
these systematic disagreements and the presence of un- 
usual HFS of the transitions to the s?Y.^ state levels - a 
phenomenon observed only in KRb up to now. Careful 
examination of the spectra and collection of new experi- 
mental observations proved that such a connection is not 
indicated by our data set. Therefore, we were forced 
to revise the initially established vibrational numbering, 
since our experience from the a^E"*" state in NaRb had 
shown that a wrong vibrational numbering could lead 
to similar systematic disagreements for high rotational 
quantum numbers. We tried several hypothesis of shift- 
ing v" by up to 2 vibrational quanta and finally for one 
of them we got an excellent description of the whole data 
set within the experimental uncertainty. 

Here we want to direct special attention to this prob- 
lem, because it is common practice to establish vibra- 
tional numbering by using the mass relations between the 
isotopes. Our experience in NaRb and KRb a'^S"'" states 
demonstrates now, that given a limited set of levels with 
assinged rotational and vibrational quantum numbers, it 
is still possible to achieve an almost satisfactory quality 
of the fit assuming an incorrect vibrational assignment. 
In both molecules, discussed as examples here, we were 
able to construct PECs that fit the experimental data 
up to N"=70, but turned out to be systematically in- 
adequate when going to higher N". It is important to 
note also that in both cases the lowest vibrational levels 
{v" = — 2 for KRb) were not observed experimentally 
thus missing in the data field. 

Finally, the analysis of our experimental data led to 
the assignment of 2400 transition frequencies to about 
1100 rovibrational levels in the a'^E"'" state in KRb. In 
Fig.mthe range of the observed triplet state energy levels 
is presented. The full list of the transition frequencies is 
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FIG. 4: (color online) The range of vibrational and rotational 
quantum numbers of the energy levels of the a^E^ state, ob- 
served in the present study. 



given in Table I of the supplementary material. 
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FIG. 5: (color online) A portion of progressions from the same 
upper state level (v', J'=51), but excited from different lower 
X^E"*" state levels: (v"=9, J"=50) for the upper trace and 
(v"=8, J"=52) for the lower trace. 



B. The X^E+ state 

After the assignment of the levels of the a'^I]+ state 
and the derivation of a preliminary potential for it, we 
made first coupled channel calculations including the two 
states correlated to the 4s-|-5s atomic asymptote and the 
hyperfine interaction The purpose was to determine 
the region where the X^E"*" and a^I]+ states could be 
treated as uncoupled, and to use this region for fitting 
the PEC for the t^Y^ state. For the XiE+ we took the 
potential given in Ref. [13], extended it by the dispersion 
coefficients from Refs. |9|, ilOj. These calculations indi- 
cated that the mixing between the two states is expected 
to be strong around = 86 and exactly there we found 
the largest disagreement between the results of the sin- 
gle channel fit and the experimental observation. This 
disagreement, however could not be fully explained even 
after the coupling was taken into account. Therefore, we 
decided to reinvestigate experimentally the highest lev- 
els of the ground state in KRb, trying to collect higher 
vibrational levels with a wider range of J". 

Due to better signal-to-noise ratio in our experiment 
we were able to extend the highest observed vibrational 
quantum number to v"=88. The corresponding levels 
were observed for J" = 51 — 55. The assignment of 
the observed transitions was based on the analysis of 
Ref. 13]. As already mentioned in Section |lT] the weak 
KRb progressions close to the asymptote were overlapped 
by stronger K2 transitions. Apparently, this led to an in- 
correct determination of the transition frequency to the 
(v"=86, J"=52) level in Ref. [H by about 0.02 cm"!. 
In order to increase the confidence in our assignment 
and avoid similar situations we excited the same up- 
per level from different X^S"*" levels. Thus, KRb lines 



of interest should be the same, but overlapping KRb or 
K2 lines should be shifted or disappear. As an exam- 
ple, we show in Fig. [5] progressions starting from (v', 
J'=51), excited from (v"=8, J"=50) for the upper trace 
and (v"=9, J"=52) for the lower trace. The same pro- 
gression is shown also in Fig. 1 of Ref. [3], but there 
excited from (v"=9, J"=52). By a simple inspection it is 
immediately clear what lines appear in both traces and 
are thus the desired KRb lines. Overall, the present ex- 
periment brings 1200 additional lines to the data set of 
the XiS+ state. 

In order to make a reliable connection in the energy 
scale between the X^S+ and the a^S+ states we added 
to the ground state data set also levels of progressions 
(about 600 transitions) that were recorded simultane- 
ously with progressions to the triplet state and originat- 
ing both from a common upper state level. In this way 
we add in total about 1800 new transition frequencies to 
the primary data set obtained in Ref. Ts]. 



IV. POTENTIALS 

The full description of the experimental data needs 
a model including adiabatic PECs of the two electronic 
states and parameters that describe the coupling between 
them due to the hyperfine interaction. Similarly to the 
approach adopted in Refs. fT, ^3] we first constructed the 
potentials using single channel fits and then we performed 
coupled channels calculations and determined the shifts 
due to the interaction between the states. As a next step, 
we corrected the experimental data with the calculated 
corrections, obtaining in this way "adiabatic" data sets, 
i.e. data sets that would have been observed if the states 
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would not couple. Afterwards the whole procedure was 
repeated until the results of the coupled channels calcu- 
lations agreed with the experimental observations. 

For the potential construction we first used the point- 
wise representation. In the initial stages of fitting it has 
the advantage to be very flexible, achieving a convergence 
in only few iterations. In this way the potential curve was 
refined following the growth of the experimental data set, 
several hypothesis of the vibrational numbering in the 
a'^S"'" state were tested and also some wrong line assign- 
ments were identified. For the coupled channels fit the 
pointwise potentials were transformed into analytic form 
[25i] with which the number of free parameters is reduced 
by almost a factor of two. 

We split the representation of the potentials into three 
regions: the repulsive wall (R<Ri„„), the asymptotic re- 
gion (R>Ro„f), and the intermediate region in between. 
The analytic form of each potential in the intermediate 
range is described by a finite power expansion with a 
nonlinear variable function ^ of intcrnuclear separation 
R: 



m 



R — Rm 

R + bRr, 



\JmiR)=J2a,aRy 



(1) 



(2) 



1=0 



where {a;} are fitting parameters and b and Rm are cho- 
sen during the transform process from the spline repre- 
sentation to the analytic one, Rm is close to the value of 
the equilibrium separation. The potential is continuously 
extrapolated for R < Rm„ with: 



UsR(i?) =A + B/R"" 



(3) 



by adjusting the A and B parameters. 

For large intcrnuclear distances (R> Ro«t ) we adopted 
the standard long range form of molecular potentials: 



UMR) = C/co - Cg/R"^ - C^/R^ - Cio/R'° ± i^exch (4) 
where the exchange contribution is given by 



Scxch = A.^R-' exp(~/3i?) 



(5) 



and Uoo is the energy of the atomic asymptote (excluding 
the hyperfine energies) with respect to the minimum of 
the X^I]+ state. It coincides with the dissociation energy 
of this state, D^. The exchange energy is repulsive for 
the triplet state (plus sign in ^) and attractive for the 
singlet state (minus sign). All parameters in Eqs. (|4I5[) 
are common for the X^S+ and the a'^E'^ states. 

As a first guess we used the theoretical estimates 
of long range parameters published by Derevianko et 



al. and Porsev et al. [i, [3 for the and the 

triplet a'^S+ state. For the exchange energy amplitude 
we adopted the estimation from theoretical potentials 
in Ref. The continuation at Rout was obtained 

smoothly by variation of Uoo or ao , the latter is preferred 
because the dissociation asymptote is taken as common 
energy reference for both states and set to zero here. 

For the coupled channels calculations the hyperfine in- 
teraction is taken into account as caused by the Fermi 
contact interaction. The effective Hamiltonian for the 
electronic spin S and the nuclear spin I defines the inter- 
action parameters: 



H — Ajih ■ Sjih ■ lub + Ak ■ Sk ■ Ik 



(6) 



From the experimental spectra the Ak and A^b con- 
stants could be determined applying Hund's case bps as 
basis. Generally, the molecular coupling parameters cor- 
respond to the atomic constants A. But for atoms form- 
ing a diatomic molecule the dependence of these param- 
eters on the bond length and thus on the vibrational 
and rotational quantum numbers cannot be excluded (see 
Ref. [1^). Within the resolution of our experiment we 
recorded partially resolved hyperfine structure of the lev- 
els of the a^E"*" state and since our data covered a broad 
range of such levels we checked the dependence of Ak, 
-'^Rb ^rid A|j^ on vibrational and rotational quantum 
numbers. As obtained in other cases of mixed alkalis 
the observed hyperfine splitting is described within ex- 
perimental accuracy by using the atomic parameters as 
compiled by Arimondo [2l|. Then, each special case of 
strong coupling between the singlet and the triplet levels 
is well described by the coupled channels calculations. 

With the definitions by equations [T] to [HI the itera- 
tion process of potential fitting and calculation of the 
corrections by singlet-triplet coupling was performed on 
all spectroscopic observations of the X-'^E"'" state by [l^ 
and from the present work and of the a'^E+state. It is 
especially important to note that we obtained a large 
set of fluorescence progressions of X-'^E+and a'^E+with 
a common excited level, which fixes the relative energy 
positions of singlet and triplet level schemes precisely 
and gives information on the exchange energy from lev- 
els close to the dissociation asymptote. In total, 12502 
transition frequencies were contained in the fit and 68 
free parameters were varied, which gave a dimensionless 
standard deviation a — 0.68 for both states together. 
Tables HI and [TTl present these parameters, but the values 
actually given there are from the final step of fitting de- 
scribed in the following section including data from cold 
collision studies of other authors. For the convenience of 
the reader also the position of the potential minimum in 
energy (Te) and intcrnuclear separation (Re) is given in 
both tables. 
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TABLE I: Parameters of the analytic representation of the 
X^E+state potential. The energy reference is the dissociation 
asymptote. Parameters with * are set for continuous extrap- 
olation of the potential. 



R < Jii„„ = 3.000 A 


A* 


-0.731001924x10"* 


cm 


-1 


B* 


0.590300435x10^ cm-^J 




i?inn <R< Rout = 11.000 A 


b 




-0.39 




4.06818602 


A 


ao 


-4217.814754 


cm 


-1 


ai 


0.5530165860346669x101 


cm 


-1 


02 


0. 1402690940409950 x 10^ 


cm 


-1 


as 


0. 1048318435177453 x 10^ 


cm 


-1 


a4 


-0.4087971898030580 x 10* 


cm 


-1 


0-5 


-0.1817659880541463x10^ 


cm 


-1 


ae 


-0.2739889053154796 x 10^ 


cm 


-1 


ar 


-0.3930730951201573 x 10^ 


cm 


-1 


as 


0.1755013192862486x10'' 


cm 


-1 


ag 


0.4366867152340260 x 10** 


cm 


-1 


OlO 


-0.6031147948206772X lO'' 


cm 


-1 


ail 


-0.7948028201366105x10'' 


cm 


-1 


0.12 


0. 1 172293247025888 x 10^ 


cm 


-1 


0,13 


0.920945525655671 1x10* 


cm 


-1 


a 14 


-0. 1533597323741656 x lO^" 


cm 


-1 


115 


-0.6982526182585652 x 10^ 


cm 


-1 


aie 


0. 1418207140035536 x lO" 


cm 


-1 


air 


0.3514113318792982 x lO^o 


cm' 


-1 


ai8 


-0.9529625435553510 X lO" 


cm 


-1 


aig 


-0.1096758619400216x10" 


cm 


~i 


020 


0.4743826569413903 x lO^^ 


cm 


-1 


021 


0. 1477139521952726 x lO" 


cm 


-1 


122 


-0.1768149926253122 xlQi-'^ 


cm 


-1 


02a 


0.3496665669778944X lO" 


cm 


-1 


024, 


0.4947896972615874X lO^^ 


cm' 


-1 


025 


-0.2373506540656795x10^2 


cm 


-1 


126 


-0.1033913349509970x10" 


cm 


-1 


027 


0.6234843170080333x10*2 


cm 


-1 


028 


0. 1588626857834081 x 10** 


cm 


-1 


029 


-0.9740574795906707x10*2 


cm 


-1 


030 


-0.1741861156434701x10*"' 


cm 


-1 


031 


0.9410306699092036 x 10*^ 


cm 


-1 


032 


0.1289614354386887x10** 


cm 


-1 


033 


-0.5220986485346664x10*2 


cm 


-1 


034 


-0.5776754651606419 x 10*^ 


cm 


-1 


035 


0. 1278431437418345 x 10*2 


cm 


-1 


136 


0.1182411100178171x10*^ 


cm 


-1 


-Rout *^ R 


Uoo 


0.0 


cm 


-1 


Ce 


0.2072097x10* cm-*7 





Cs 0.6509487 xlO^ cm"* A* 

Cio 0.2575245x10** cm-*A*o 

0.1469387x10^ cm-*A-T 
7 5.25669 
/3 2.11445 A-* 

Derived constants: 
equilibrium distance: Rf= 4.06770(5) A 



TABLE IL Parameters of the analytic representation of the 
a'^E^ state potential. The energy reference is the dissociation 
asymptote. Parameters with * are set for continuous extrap- 
olation of the potential. 



R < = 4.956 A 


A* 




-0.102755858x10"* cm"* 


B* 




0.595720017x10^ cm-*A* 




Rinn 


<R< Rout = 11.000 A 


b 




-0.41 


Rm 




5.90210000 A 


ao 




-249.030467 cm"* 


ai 




-0.8249570532852936 cm"* 


02 




0.1696361682907917x10"* cm"* 


as 




0.2195618553936388x10* cm"* 


04 




-0. 1830331320330776 xlO^ cm"* 


as 




-0.3577610246567447x10"^ cm"* 


ae 




0.7605943179592050x10"^ cm"* 


ar 




0.2321948211354688x10* cm"* 


as 




-0.2969108359066523x10* cm"* 


ag 




-0.9066958465430735x10^ cm"* 


aio 




0.1193993210361681x10*° cm"* 


ail 




0.2271299134652605x10** cm"* 


ai2 




-0.3822195877987927x10** cm"* 


0-13 




-0.3733663448209457x10*2 cm"* 


Ol4 




0.8258881465522462x10*2 cm"* 


015 




0.3983625994810529 xlO*-'^ cm"* 


016 




-0.1170134952322679x10*"* cm"* 


017 




-0.2572921814049761x10** cm"* 


018 




0.1083172601973219x10*5 cm"* 


019 




0.7316944602806914x10** cm"* 


020 




-0.6401304321393455x10*5 cm"* 


021 




0.2090674461383090x10*5 cm"* 


022 




0.2236346785424132x10*" cm"* 


023 




-0.2555345293959939x10*'^ cm"* 


024 




-0.3556041967731776x10*8 cm"* 


025 




0.8321268009266810x10*8 cm"* 


02(i 




-0.1364612053210692x10*8 cm"* 


027 




-0.9030566262973000x10*8 cm"* 


028 




0.8950177388640356x10*8 cm"* 


029 




-0.2723988913077922x10*8 cm"* 


Rout < R 






0.0 cm"* 


Ce 




0.2072097x10* cm-*A8 


Cs 




0.6509487x10^ cm-*A* 


Cio 




0.2575245x10** cm-*A*o 


Aex 




-0.1469387x105 cm-*A— >" 


7 




5.25669 


/3 




2.11445 A-* 


Derived constants: 


equilibrium distance: 


R^= 5.9029(1) A 


electronic term 


energy: 


r°= -249.031(10) cm-* 
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V. COLD COLLISIONS 

As it was already mentioned in the introduction, Fes- 
hbach resonances for ^°K*^Rb were observed by several 
groups 0, @| • In this section we apply the derived po- 
tentials to calculate these Feshbach resonances and com- 
pare them with observations. It is important to note 
that the spectroscopic observations of the present study 
were done on different isotopomers, thus application of 
the derived potentials for the Feshbach resonances as- 
sumes the validity of the Born-Oppenheimer approxima- 
tion within the desired accuracy. We use for the compar- 
ison the s-wave resonances observed by Rcf . [5] , because 
these data represent the largest set and come from the 
same lab, thus they should not show possible internal in- 
consistency like calibration errors of the magnetic field. 
All calculated resonances were found within 50 Gauss 
from the observed ones, this deviation is not surprisingly 
large because these calculations are extrapolations out of 
the range of spectroscopic observations; the largest outer 
turning point of observed levels is at 12.6 A for the sin- 
glet state and at 14.8 A for the triplet state being 6.4 
cm~^and 2.3 cm^^below the dissociation asymptote, re- 
spectively. This indicates, that the spectroscopic obser- 
vation ends within the changeover to the long range be- 
havior and the quality of the extrapolation follows more 
from the reliability of the theoretical estimates of the 
dispersion coefficients [lol | than from our experimental 
data. 

Thus the measured Feshbach resonances contain infor- 
mation on the long range function. To include them in 
the total potential determination we set up a non-linear 
least squares fitting code for the dispersion parameters, 
which contains the scattering calculations as subroutine 
and searches for the maximum of the elastic rate coeffi- 
cient which corresponds to the observed Feshbach reso- 
nance. The resonances have been observed by trap loss of 
potassium, thus by an inelastic process, which will need 
three body collisions. As very often assumed, the three 
body cross section should be enhanced by the condition 
of a two body Feshbach resonance. The calculations were 
done for a temperature of IfiK and use the accuracy of 
the magnetic field of 0.2 G, as stated in Ref. [il, as uncer- 
tainty for each individual observation. In a first step the 
resonances were fitted by varying the dispersion terms 
Cg and Cg, only. With the new total potentials asymp- 
totic levels of pure singlet (v — 97, 98, 99) and triplet 
(v=:29, 30 31) states were calculated and used as input 
for a new fit of the single channel calculations. Error 
limits for these levels were estimated from the fit of the 
resonances to be about 100 kHz. In the new fit of the 
inner part of the potentials the coefficient Cio was cho- 
sen as an additional free parameter. With this result a 
new iteration was started. Only two such iterations, res- 
onance fit followed by single channel fit, were needed for 
convergence. The final results are given in table [T] for 
XiE+ and in TableHHfor a.^E+. 

The ten s-wave Feshbach resonances could be fitted to 



a standard deviation of cr = 1.37 and show on average 
a systematic deviation of -1-0.05 G, i.e. the observations 
are at larger magnetic fields than the results from the fit. 
This fit of the resonances is slightly better than the one 
published in j5j: cr — 1.66 and mean deviation -1-0.04 G, 
and has the additional advantage, that all spectroscopic 
observations are equally well described as stated in the 
previous section. All long range parameters are close to 
the theoretical predictions [§,[l3| , which indicates that we 
have obtained a physically consistent model for the two 
electronic states correlated to the atomic ground states 4s 
K -f 5s Rb. It will be applied now to calculate an overview 
of expected collision properties at low kinetic energy, the 
upper limit of which is not easy to estimate because of 
the artificial repulsive branches of the potentials (see eq. 

El). 

Ferlaino et al. Q give a list of singlet and triplet 
scattering lengths for all possible isotope combinations 
of KRb with natural abundance and calculate the scat- 
tering length of the lowest hyperfine input channel of 
KRb. For the application of mass scaling they used full 
potentials from other sources which give for ^°K^^Rb as 
number of supported vibrational levels 98 for singlet and 
32 for triplet. The potentials derived in the present work 
support also 32 for triplet, however 100 for singlet. Such 
possible variation was taken into account by the error 
limits stated in |4j]. Table Hill compares our results with 
those from Ferlaino et al. The two columns, labeled as 
o-iowest, show scattering lengths of the lowest Zeeman 
level in each case of the isotope pair, e.g. for ^^K -|- 
^^Rb fA'=l, mK=l and fi^h=2, mRb^2. 

Certainly, for ''''K^^Rb the results closely agree, but 
for the other isotopomers significant differences, espe- 
cially for the singlet case appear, which reflects the dif- 
ferent number of supported vibrational levels within the 
respected potentials of this state. 

In the experiments by Ferlaino et al Q and Ospelkaus 
et al '6'| also p wave resonances and s wave resonances 
coupled to d levels are reported. We checked if these are 
also well described with the present potentials. Espe- 
cially in @ the splitting between mi = and to; = ±1 
of a p wave is reported and it gives the opportunity to 
check if the spin-spin or second order spin-orbit coupling 
for the triplet state needs to be adjusted which both to- 
gether can be described by the effective operator ^t\ : 

Hss^lxi3S^,-S'). (7) 

A is a function of internuclear separation R and typi- 
cally is derived from spectroscopic data as an expectation 
value of a specific vibrational level. Here we use simple 
functional forms to incorporate the pure dipole-dipole 
contribution for the spin-spin interaction as a function 
and the spin-orbit part as an exponential function 
refiecting effectively an overlap integral of the electronic 
distribution of the two atoms: 
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TABLE III: Comparison of scattering lengths (in Bohr radius, [ao — 0.5292 x 10 ^"m). 



isotope 




glet 


^triplet 




^lowest 






present 


\5] present 




present 


39/85 


26.5(9) 


33.4 


63.0(5) 


63.9 


56.6(4) 


58.0 


39/87 


824l?« 


1868 


35.9(7) 


35.90 


27.9(9) 


28.36 


40/85 


64.5(6) 


65.8 


-28.4(16) - 


28.55 


-25.3(16) 


-21.12 


40/87 


-111(5) 


-111.5 


-215(10) - 


215.6 


-185(7) 


-183.1 


41/85 


106.0(8) 


103.1 


348(10) 


349.8 


283(6) 


283.6 


41/87 


14.0(11) 


7.06 


163.7(16) 


164.4 




647.8 



As5 = -ja^ {^-^ + aso exp {-b{R - Rso))j (8) 

All quantities are given in atomic units, a is the fine 
structure constant and aso, b, and Rso are the model 
parameters for the second order spin-orbit contribution. 
Since the data in hand cannot be sensitive to the actual 
function, we selected values for b and i?so, which are 
comparable to the example of Rb2 28] (b= 0.7196 a^^ 
and Rso= 7.5 ao) and fitted the parameter aso which 
gives aso= -0.013 ag^ (ao = 0.5292 x IQ-^" m). The 
result shows that the first part in eq. [5] dominates in the 
internuclear separation interval R > 12ao important for 
the bound levels of the triplet state which determine the 
Feshbach resonances. 

Taking this result into account, the four resonances 
reported in d, Q involving / — 1 and 2 were calculated. 
The I = 1 resonances were found within the experimental 
accuracy, and it is predicted that the resonance at 456G is 
split for m;=0, ± 1 by about 0.1 G where the component 
m;=0 is higher in field than the others. This is opposite 
to the case of the observed resonance at 515 G. The I = 2 
resonances were found above the observed magnetic field 
but not more than 0.7 G, the reason for this difference 
is presently unclear. These resonances are very sharp, 
the calculated rate constants of the corresponding elastic 
collisions have half widths of less than 0.001 G, which 
depends on the magnitude of the above mentioned SS 
coupling. No value was given in [5] for their calculations 
of the widths. 



VI. CONCLUSIONS 

An extensive dataset on the a^E+ state of KRb was 
collected covering energy levels in a wide range of vi- 
brational and rotational quantum numbers in ■^^K^^Rb 
and ^^K^^Rb. The classical turning point of the last ob- 
served level iv' = 25, J' = 26) hes at 14.6 A on the fitted 
rotationless potential energy curve. The long range pa- 
rameters Cg, Cs, and Cio deviate only slightly from the 
theoretical estimates by ref. [1, , but we should remind 
the reader, that the present data set does not allow to 
determine independently all three parameters. Thus due 



to correlations between those of the data analysis, there 
is no contradiction to those of the theory. But for mod- 
eling of cold collisions one should use the values reported 
in Tables HI and [TTl of the present work. 

The complete determination of the potential curves al- 
lows for reliable extrapolation to the dissociation limit 
and thus their depths. Tables U and [ll] give 4217.815(10) 
cm~^ and 249.031(10) cm^^ for the respective dissocia- 
tion energies. Recently, Wang et al. [l^ reported about 
an independent determination of the dissociation energy 
of the singlet state using a photoassociation-deplction 
process to obtain binding energies of loosely bound vi- 
bration levels, namely w = 87 and 89. The dissocia- 
tion energy was derived by adding level energies with re- 
spect to the minimum of the potential curve as reported 
in Ref. by a near-dissociation expansion (NDE) ap- 
proach. The final value 4217.822(3) cm~^ agrees with 
ours within the given error limits. But for estimating 
the error limit Wang et al. assume that the accuracy of 
level energy w = 87 is equal to the standard deviation 
of the NDE fit, 0.0015 cm~-^. Looking into the actual 
data set from Ref. [l3| one finds that for w = 87 only two 
high rotational states with J = 50 and 52 were observed. 
Thus a fairly far extrapolation to the levels with J = 
is done, and one should certainly include the correlation 
of the NDE parameters for the error estimation. Addi- 
tionally, we would remind the reader, that we corrected a 
misidentification of v=86 levels, due to overlapping lines, 
and this would effect the NDE analysis. Therefore, the 
good agreement between our present dissociation energy 
and that reported by Wang et al. could be fortunate. 
As a second check we calculated the binding energies for 
f = 87 and 89 from our results and found them larger 
by 0.026 cm~^ and 0.053 cm~^, respectively. For u = 87 
the deviation is about 12 times the error limit reported 
by Wang et al. We tried to incorporate the measured 
binding energies in a new fit of the potentials including 
our spectroscopic results and the Feshbach resonances 
from Ref. but no convincing result was obtained up 
to now. The reason for the discrepancy is unclear for us, 
especially because we extended the dataset from Ref. [ll] 
just around w = 87 with more J levels and by measure- 
ments with u = 88. Thus the energy interval where the 
w = 87 is located is well covered by our measurements 
which have an accuracy of typically 0.003 cm~^. 
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By new spectroscopic data the potential curve of 
X^S+could be refined compared to the work by Ref. 
[T^ . and we were able to link it to the potential curve 
of a^E+with respect to the asymptotic function and to 
the relative energy position of both functions. This was 
the prerequisite to start the inclusion of the observed 
Feshbach resonances Q into the modeling, which finally 
leads to a complete system for predicting ultracold col- 
lision properties. Such prediction can be extended to 
other isotopomers, as presented partly in Table Hill But 
it would be clearly desirable to get precise Feshbach res- 
onances also for different isotopomers in order to check 
the degree of validity of the Born-Oppcnheimer approxi- 
mation. 

Despite the combination of the spectroscopic and Fesh- 
bach data there is a significant gap in the energy interval 
below the atomic asymptote (about 3 cm^^) in which 
bound levels of low rotational quantum numbers, prefer- 
able J or N =0, could be observed. Thus the predic- 
tion in this region relies on the quality of the long range 
parameters and measurements e.g. by two-color photo 
association would be valuable. To reach such levels by 
Feshbach spectroscopy is not promising because too high 
magnetic fields would be needed. 

KRb is very often discussed for applications in stud- 



ies of dipolar cold gases or even quantum computing. To 
transfer Feshbach molecules to deeply bound levels which 
are stable and would have large electric dipole moments, 
spectroscopic data on excited states will be of importance 
to predict efhcient processes. Such data are partly exist- 
ing from other authors and arc in hand also from our 
own spectroscopic work, which will be analyzed joining 
the available data. Of particular interest in KRb is the 
HFS in the B^II state which is supposed to cause the 
"irregular" hyperfine splitting the the transitions to the 
a^I]+ state. 
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